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ABSTRACT 

More than 50% of Active Galactic Nuclei (AGNs) are suspected to be red and affected by dust- 
obscuration. Meanwhile, popular spectral diagnostics of AGNs are based on optical or ultraviolet light, 
making the dust obscuration as a primary concern for understanding the general nature of AGNs and 
supermassive black holes residing in them. To provide with a method of investigating properties 
of the dusty AGNs, we derive new black hole (BH) mass estimators based on velocity widths and 
luminosities of Near Infrared (NIR) hydrogen emission lines such as Pa and P/?, and also investigate 
the line ratios of these Hydrogen lines. To derive the BH mass (Mbh) estimators, we used a sample 
of 37 unobscured Type-1 AGNs with a Mbh range of lO^'^-lO^ '' Mq, where Mbh come from either 
reverberation mapping method or single-epoch measurement method using Balmer lines. Our work 
shows that Mbh can be estimated from the Paschen line luminosities and the velocity widths to the 
accuracy of 0.18 - 0.24 dex (rms scatter). We also show that the mean line ratios of the Paschen lines 



and the Balmer lines are 



9.00 



M 

' Pc 



2.70, which are consistent with a Case B recombination 



under a typical AGN broad line region environment. These ratios can be used as reference points 
when estimating the amount of dust extinction over the broad line region (BLR) for red AGNs. We 
expect the future application of the new BH mass estimators on red, dusty AGNs to provide a fresh 
view of obscured AGNs. 

Subject headings: black hole, AGN, dusty AGN, NIR, Paschen line 



I. INTRODUCTION 

Supermassive black holes (SMBHs), found at centers 
of massive spheroids and active galaxies, are considered 
to play an important role in the formation and the evo- 
lution of galaxies. It is suggested that SMBHs reg- 
ulate the star formation activities of galaxies through 
their enormous energy output during active phase, pro- 
viding a feedback mechanism to reconcile the observed 
trend of downsizing of the galaxy evolution with hier- 
archical galaxy formation models (lJuneau et al.l 120051 : 
iBundv et al.l l2006':'Sch awinski et al.ll2006( ). However, the 
coevolution of SMBHs and their host galaxies stand as an 
unsolved astrophysical problem, as such a process needs 
to eventually lead to a rather unexpected tight corre- 
lation between host galaxy properties and SMBH mass 
todav (iFerrarese fc MerrittI l2000t iGebhardt et"all 120001: 
iTremaine et al.ll2002f ). 

One of the great challenges toward the understand- 
ing of the evolutionary sequence of AGNs and the 
co-evolution of SMBHs and their host galaxies comes 
from the so-called red or dusty AGNs which are be- 
lieved to occupy more than 50% of the AGN population 
(iCoinastri et al.l 120011: iTozzi et al.l 120061: iPolletta et all 
|2008[ ). Red, dustv AGNs are expected, if the initial phase 
of AGN activity starts in a dust-enshrouded environment 
such as within massive starburst regions of luminous in- 
frared galaxies. In such a scenario, AGNs become more 
visible after sweeping away cold gas and dust that are 
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necessary to s ustain the massive star form ation in their 
host galaxies (jHopkins fc Hernauistl[2006l) . On the other 
hand, dust obscuration can arise in a different way in the 
unified model of AGN when accretion disks and broad 
line regions (BLRs) around SMBHs are viewed through 
dust torus. In any case, red, dusty AGNs can shed light 
on the properties and the evolution of the AGN popula- 
tion in general. 

The phenomcnological definition of "red, dusty AGNs" 
is rather broad. It can include AGNs selected in a variety 
of ways such as those selected by very red colors in opti- 
cal through NIR and the radio detection (e.g., R — K>5 
mag and J — K > 1.3 mag of the sample in Urrutia et al. 
2009; also see Cutri et al. 2001), red MIR colors (Lacy 
et al. 2004; Lee et al. 2008), and hard X-ray detections 
(e.g., Polleta et al. 2007). These different kinds of AGNs 
have a common characteristic where their SEDs are red, 
due to the obscuration of their light by the foreground 
gas and dust. Hence, they are considered to be the in- 
termediate population between the dust-enshrouded star 
forming galaxies and the unobscured AGNs. The dust- 
obscuration does not necessarily exclude Type-1 AGNs 
(AGNs with broad emission lines; e.g., Alonso-Herrero 
et al. 2006), as red, dusty AGNs with broad emission 
lines are found quite often. Of 23 X-ray QSOs with IR 
power-law SED and spectroscopic redshift identification, 
14 are classified as broad-line AGNs (Szokoly et al. 2004; 
Alonso-Herrero et al. 2006). More than 50% of the red 
AGNs with the radio and red optical/NIR color selection 
are found to be type-1 AGN which can have the redden- 
ing parameter of 2 or more (Uruttia et al. 2009; Glikman 
et al. 2007). The existence of broad line AGNs among 
red, dusty AGNs opens a possibility of studying the prop- 
erties of the dust-obscured quasars in more detail using 
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traditional type-1 AGN diagnostics. 

However, even the measurement of the most basic AGN 
parameters such as Mbh and the Eddington ratio (the 
accretion rate) is a difficuh task for the red, dusty AGNs. 
In general, A/bh are derived from the optical or the ultra- 
violet (UV) part of AGN spectra for which spectral diag- 
nostics are well established. In the case of dusty AGNs, 
the dust obscures or significantly reduces the UV or the 
optical light coming from the region around SMBHs, 
making the popular AGN optical/UV spectroscopic di- 
agnostics useless. For example, popular A/bh estimators 
are based upon a virial relation between two parameters 
- the velocity width of the broad. Ha or H/3 lines, and 
the size of BL R estimated from t he continuum luminos- 
ity at 5100 A (iKaspi e t al.ll2000|) or the luminosities of 
Ha or H/3 ([Gree ne & Hol l200,'5n . If the light from a red, 
dusty AG N is extincted by a color excess oi E(B — V) 
=2 mag (jGlikman et all 120071 : lUrrutia et a"ni2009D . its 
Ha and H/3 line fluxes would be suppressed by a factor 
of 100 and 1000 respectively. One can try to estimate 
the amount of the dust extinction from a continuum flt- 
ting of the optical-UV spectrum or through the Balmer 
decrement, but such estimates are often inconsistent with 
each other and accompanied with uncertainties of order 
of 5E{B — V) ^ 0.5 mag or more which are related to 
the dispersion in the intrinsic properties of AGNs. The 
dust obscuration can arise from the Galacitic extinction 
for AGNs at low galactic latitude (e.g., Im et al. 2007; 
Lee et al. 2008). 

The problem can be substantially alleviated if we can 
use NIR lines instead of optical or UV lines, since NIR 
Hydrogen lines such as the Pa and P/3 are much less af- 
fected by the dust extinction than the UV/optical light. 
For the red, dusty AGN with the color excess of E{B — V) 
= 2 mag, the line fluxes of Pa and P/3 are suppressed 
by a factor of only 2.3 and 4.7 respectively. This is 
a signiflcant improvement over the optical lines. Since 
AfsH oc L°-^, where L is a luminosity of the continuum 
or an emission line, the suppression in the Paschen line 
luminosities introduces uncertainties in A/bh estimates 
only at the level of a factor of 2 or less, even without 
correcting for the dust extinction with the Paschen decre- 
ment. 

In this paper, we will derive the A/bh estimators based 
on the NIR Hydrogen lines of Type-1 AGNs and investi- 
gate the line ratios of the Paschen lines, keeping in mind 
future applications of such relations to studies of dusty, 
red AGNs with broad emission lines. 

2. DATA 

2.1. The Sample 

In order to construct mass estimators based on the 
NIR Hydrogen lines, we used two samples of low redshift 
Type-1 AGNs with NIR spectra at 0.8-4.1 jim. One 
is from Landt et al. (2008; hereafter L08) who stud- 
ied 23 well-known Type-1 AGNs in the local universe. 
Among these, we use 16 and 21 AGNs that have line 
flux and width measurements in Pa or P/3 lines, respec- 
tively. Note that we excluded four objects, Mrk 590, 
NGC 5548, Ark 564, and NGC 7469 in all or some of our 
analysis even though Pa and P/3 line flux and FWHM 
measurements were available. For Ark 564 and NGC 
7469, the FWHM values of Pa and P/3 lines listed in 



L08 are too discrepant from each other differing by more 
than a factor of 1.5. This suggests that one of the two 
measurements is erroneous. Therefore, these two objects 
are excluded from all of the analysis performed below. 
For Mrk 590 and NGC 5548, the Ha and the H/3 widths 
listed in LOS differ by more than a factor of 1.5, but 
the Paschen line widths are consistent with each other. 
These two objects are excluded in the analysis of the cor- 
relation between the Paschen and the Balmer lines, but 
their reverberation-mapping derived A/bh are retained 
for the derivation of the black hole estimators (of the 
method 2 and the method 3 in Section 3.2). 

Another sample comes from Glikman et al. (2006; 
hereafter G06). The G06 sample is made of 26 Type-1 
AGNs that were selected from a cross- match betwee n the 
SDSS-DRl quasar catalog (Schneid er et al.l I2003D and 
the 2 micron AU Sky Survey (2MASS; Skrutskie et al. 
2006), with the criteria of K < 14.5 mag, redshift z < 
0.5 and the absolute magnitude in z-band. Mi < -23 mag. 
Here, the absolute magnitude limit is imposed to mini- 
mize contamination of the quasar light due to the host 
galaxy. The G06 sample provides 11 AGNs with Pa lines 
and 14 AGNs with P/3 lines showing up in their NIR 
spectra. 

In all, we use 27 and 35 Type-1 AGNs for the Pa and 
the P/3 line analysis, respectively. Table [1] summarizes 
properties of these AGNs. 

2.2. Black Hole Mass and Analysis of Optical Spectra 

For A/bh of our sample AGNs, we use the following 
values. 

First, we use A/bh derived with the reverberation map- 
ping method from literatures if available. This applies to 
10 AGNs with Pa and 11 AGNs data in L08 for which 
A/bh are taken from Vestergaard & Peterson (2006). 

Second, if A/bh from the reverberation mapping 
method are not available, we derived A/bh from the op- 
tical spectral information using single epoch A/bh esti- 
mators. For the L08 sample without the reverberation 
mapping derived A/bh, we used the H/3 luminosity and 
FWHM hsted in the Table 5 of L08 to derive their A/bh- 
For this, we adopted a AIbh estimator in Greene & Ho 
(2005), after adjusting their relation to have t he virial 
factor of 5.5 (|Onken etal .112001 1 Woo eFalllMoh by mul- 
tiplying their relation by a factor of 1.8. We also made 
a small correction to the L08 FWHM and luminosity 
values, since the way how Greene & Ho (2005) derived 
FWHMs and luminosities is different from what LOS 
did - Greene & Ho (2005) adopted a multiple Gaussian 
component fitting to the broad line while LOS adopted 
a single component fitting of the broad line after re- 
moving narrow line components. We find that the cor- 
rection factor to be FWHM,„„/ti/FWHMs„g;e=0.9 and 
Lmuiu/Lsingie = i-OS On average through fitting of the 
SDSS spectra of the G06 sample with the two methods 
(See below). 

For the G06 sample, none of the objects have A/bh 
derived from the reverberation mapping method. Some 
of the G06 sample have their H/3 line luminosities and 
FWHMs listed in Shen et al. (2008), but after check- 
ing FWHM values in Shen et al. (2008) and the optical 
spectra, we find cases where FWHMs appear to be over- 
estimated. Therefore, we used our own spectral fltting 
routine (Kim et al. 2006) to obtain the H/3 FWHMs and 
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the H/3 luminosities from the SDSS spectra. The fitting 
model includes components for a power law continuum, 
the host galaxy light, broad Fe II multiplet emissions, 
a Balmer continuum, and multiple Gaussian profiles for 
broad and narrow emission lines. The FWHM and the 
line luminosities are derived from the best-fit model pro- 
files of the broad line component which can be consisted 
of a sum of multiple Gaussian components. 

From the derived luminosity and FWHM values, we 
estimated Mbh by using an updated version of the H/3- 
based virial mass estimator (Greene & Ho 2005). We 
updated Equation (7) in Greene & Ho (2005), by adopt- 
ing a newer i?BLR --^5ioo relation, i.e.. Equation (2) of 
(|Bentz et al.ll2009D . The updated relation IS given be- 
low. 



M 



_ q6.88±0.57/ 

Mq " ^1042 erg s 



0.46±0.05 



FWHMh^ 2 
103 km s-i-^ ■ 



Again, a proper adjustment was made for the virial 
factor in order to make the relation consistent with the 
virial factor in Vestergaard & Peterson (2006). The de- 
rived BH masses are listed in Table 1. 

2.3. Analysis of NIR Spectra and Line Information 

The G06 data were taken mostly with the SpeX in- 
strument on the IRTF with the spectral resolution of 
^ ~ 1200 at 0.8 - 2.5 fim, and ^ ~ 1500 at 1.9 - 4.1 
/im. G06 provide with the flux-calibrated, reduced NIR 
spectra, and we used the following procedure to measure 
the Paschen line fluxes and FWHMs of the G06 sample. 

The G06 spectra were converted to the rest-frame, and 
then the continuum around the Paschen lines were de- 
termined by fltting a linear function to the continuum at 
regions near the Paschen lines. We varied the wavelength 
centers and widths of the continuum fltting regions, but 
typical values are chosen to be ^2 FWHM of the broad 
component of H/3 line away from the Paschen line center 
and the width of 0.02 /im for the continuum fltting region 
centers and widths respectively. For the Pa line, this cor- 
responds to the wavelength regions centered somewhere 
at 1.82 to 1.84 /xm and 1.90 to 1.92 //m with the width 
of 0.02 fim, while for the P/3 line the wavelength cen- 
ters are somewhere at 1.23 to 1.25 fj,m and 1.31 to 1.33 
fim. When the continuum shows data points affected by 
imperfect sky subtraction, we adjusted the continuum 
fitting regions. Also, the width of the continuum fitting 
regions are varied from 0.01 to 0.03 /im, and the cen- 
ter wavelengths are shifted by 0.01-0.02 /xm to provide 
robust measurements. Although NIR spectra of Type-1 
AGNs are known to be dominated by a black body radia- 
tion from warm dust, we find that the linear fitting of the 
continuum over the limited wavelength range provides a 
reasonable approximation of the continuum around the 
Paschen lines (see Figure 2; also see L08). 

After the continuum subtraction, we used a single 
Gaussian function to fit an emission line. The fitted pa- 
rameters are multiplied by correction factors to correct 
for systematic errors in the single component fitting (see 
next paragraph). During the fitting, we set free the cen- 
tral wavelength of each Gaussian component. Therefore, 
the free parameters of the fit are FWHM, flux, and the 
central wavelength for Gaussian component. The flt was 



performed with the MPFITEXPR procedure of the IDL. 
Figures 2 and 3 show the line-fltting results. The flts 
provide a formal error of 12% in flux and 9% in FWHM. 
We also estimated the error arising from uncertainties 
related to how the continuum is determined. This was 
done by varying the regions for the continuum fitting as 
described earlier. We find that the average uncertainty 
arising from the continuum determination is 4% (fiux), 
3% (FWHM), ahhough we find cases where the FWHM 
values change by as much as 25%. Such cases occur when 
there is an extended feature at the wing of the emission 
line which may or may not be a true feature. The total 
errors in fiux and FWHM measurement are taken as a 
square root of the quadratic sum of the formal fitting 
error and the error from the continuum determination, 
which are found to be typically 12% and 9%. 

Fitting a multiple component line with a single 
Gaussian function can produce systematic errors in the 
fitted parameters, and we tried to correct for such a bias 
by applying correction factors that are derived from 
Balmer lines. This is done by comparing the broad line 
widths and the line fluxes from multiple component 
flts of Balmer lines versus the fltted parameters from a 
single component flt of the same Balmer lines. By doing 
so, we correct potential systematic arising from not 
removing narrow line components in the fltting process 
too. The correction factors are computed for each object 
and they range over 0.74 - 1.01 with the mean of 0.91 
for the FWHM ratios (FWHM„„itj/FWHMsj„gZe) and 
0.98 - 1.17 for the line flux ratio with the mean of 1.06 
{Lrnuiti I L single)- This procedure makes an assumption 
that the line proflle shapes of the Paschen lines follow 
those of the Balmer lines. To conflrm this assumption, 
we compared the P/3 line proflles with the model-fltted 
line proflles of Ha. Only the P/3 lines with the S/N 
> 25 are shown in Figure 4. The P/3 lines are plotted 
as histograms with error bars, while the Ha lines are 
indicated with thick solid lines. The thin dotted lines 
and the thick dashed lines represent the broad and the 
narrow components of the Ha lines. The P/3 proflles 
are normalized to have the maximum value of 1 while 
the Ha are normalized to have the same total flux 
as the normalized P/3 proflles. Figure 4 demonstrates 
that the P/3 proflle shapes are consistent with those of 
the Ha proflles, except for the width and the overall 
normalization, supporting our assumption that the 
Paschen line proflles are similar to those of the Balmer 
lines. We also tried fltting the Paschen line proflles with 
double Gaussian proflles for objects that seem to have 
narrow line components (SDSSJ010226.3-003904.6, 
SDSSJ015530.0-085704.0, SDSSJ015950. 2+002340. 8, 
SDSSJ031209. 2-081013.8, SDSSJ015910.0+010514.5, 
SDSSJ235156. 1-010913.3, SDSSJ000943. 1-090839.2 
and SDSSJOllllO.0-101631.8). For H/3, we flnd that 
the double Gaussian fltting returns FWHM values 
greater than those from the single Gaussian fltting 
by a factor of 1.08 in average. However, the double 
Gaussian fltting (one for narrow line, and another for 
broad line) tends to return FWHM values greater than 
the results from the multiple component fltting. Thus, 
another correction factor is needed to convert the double 
Gaussian fltting result to the multiple Gaussian fltting 
result, and the result from this double Gaussian fltting 
test should be taken seriously. We are presenting a result 
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from this test in order to see if we find the same kind 
of the dependency of fitted parameters on the fitting 
method which we find in the Balmer fines as a way to 
support our assumption that the Paschen fine profiles 
are similar to the Balmer line profiles. 

In most case, the fitted values from the double com- 
ponent fit agree with those from the single component 
fit (without the additional correction factor), where the 
FWHM and fiux values from the double component fit- 
tings are in average 7% larger and 2% smaller than 
the single component fitting results respectively. This 
agrees well with the trend we find from our analysis 
of the Balmer lines using single, double, and multiple 
component Gaussian fittings (see above). One exception 
is SDSSJ015950.2-t-002340.8 whose FWHM changed by 
about 50%. Distinguishing the narrow and the broad 
components is difficult for this object, therefore the mean 
value between the parameters from the two different 
methods was adopted as the best-fit value, with the half 
of the difference as its error. We note that the double 
component fit improved the reduced values signifi- 
cantly (> 0.3) for only three objects, where the improve- 
ment came from fitting of broad extended wings which 
did not affect the derived fitting parameter values rather 
than through the change of the FWHM values (except 
for SDSSJ015950.2-h002340.8). 

Finally, the measured FWHMs were corrected for 
the instrumental resolution and the fitting methodol- 
ogy (multiple component fit versus single component fit), 
and the fluxes were converted to the luminosity assum- 
ing a standard ACDM cosmology of Hq—70 km sec~^ 
Mpc-i, rJ™=0.3 and flA^O.? (e.g., Im et al. 1997). The 
measured line luminosities, fluxes, and FWHMs are pre- 
sented in Table [H 

LOS present the results of their line analysis, based 
on the NIR spectra that were taken with the Spex 
spectrograph (|Ravner et al.l 120031 ) on the Infrared Tele- 
scope Facility (IRTF) at an average spectral resolution 
of 400 km s~^. For the LOS sample, we use the line 
fluxes and FWHMs derived by them and presented in 
their Table [1] after correcting FWHM values for the 
instrumental resolution. We also applied the correc- 
tion factor of FWHM,„„itj/FWHMsi„gie6road=0.9 and 
Lrnuiti/ Lstngiebroad='i--0& which corrccts for the differ- 
ence in the line-fitting methods (LOS versus Greene & 
Ho 2005) to derive the line parameters which are sim- 
ilar to the correction factors we derived for the fitting 
process of the G06 spectra. As for the accuracy of the 
LOS measurements, LOS quote a typical error of 10% of 
less. Therefore, we adopt a conservative value of 10% 
for the measurement error of the fiuxes and the FWHMs 
presented in LOS. 

We note that contamination from the host galaxy light 
to these measurements is negligible. The contamination 
of the line flux due to the host galaxy is possible, but the 
G06 data were taken with a narrow slit to minimize the 
host galaxy light to less than S%. We expect that the 
same statement holds true for the LOS sample whose data 
were taken with a narrow slit for AGNs at z < 0.1 which 
are located much closer to us than the G06 sample. 



3. RESULTS 



3.1. Empirical Relation between Balmer and Paschen 

Lines 

Before constructing mass estimators based on the 
Paschen lines, we show here that how well the proper- 
ties of the Paschen lines correlate with the Balmer lines. 
A tight correlation between the two lines would imply 
that the Paschen lines originate from the broad line re- 
gions similar to the Balmer lines, thus serving as a strong 
justification for the use of the Paschen lines as a mass es- 
timator. Good correlations between the FWHM values 
of the broad Balmer lines and the broad Paschen lines 
were shown in LOS. Here, we extend the analysis to the 
line flux ratios, and add quasars from the GOG sample to 
the LOS sample to strengthen the goodness of the FWHM 
correlation. Furthermore, we also derive equations that 
relate the properties of the Balmer and the Paschen lines. 
Figure 5 shows the correlation between FWHM values of 
the Balmer and the Paschen broad lines, while Figure 6 
shows a correlation of line luminosities of the broad lines. 
To derive the correlations between the two quantities, we 
performed a linear bisector fit using the equations below. 



log( 



FWHMy 
1000 km s- 



, , FWHMx , 
^ + ^^°g( lOOOkms-i ) 



Here, X and Y are line identifiers, and A and B are 
the correlation coefficients in the fit for FWHM, and C 
and D are the coefficients for the line luminosity ratio 
fit. The results of the fitting are summarized in Table 
2. The table also lists the rms scatter of the data points 
with respect to the best-fit lines. 

These results show that the Paschen line luminosities 
and FWHMs correlate well with those of the Balmer 
lines. The rms scatters in the line luminosity correla- 
tion are ■^O. 12-0. 14 dex against Ha, and ^0.16-0.19 dex 
against H/3. For the FWHM values, the rms scatters 
are 0.045-0.06 dex against Ha, and 0.05-0.06 dex against 
H/3. The slightly larger scatters and a notable offset in 
FWHM values of H/3 against Paschen lines suggest the 
complexities in AGN spectra around H/? line noted in 
LOS, an excess, broad component in the red part of the 
H/3 line caused by an unclear origin (e.g., Meyers & Pe- 
terson 19S5; Veron et al. 2002). We also point out that 
the line widths of Balmer lines are systematically larger 
than those of the Paschen lines, and that the trend is 
stronger as the wavelength decreases. This suggests that 
the Paschen broad lines and the Balmer broad lines orig- 
inate from a similar BLR, but with Balmer lines origi- 
nating from the inner region of the BLR than Paschen 
lines. 

Similarly, we also examine correlations between 
^^5iooA-^ and Lpo,,^. -Rblr values are derived from 
L(5100) using Equation (4) of (IGreene fc Hol[2005l) . Fig- 
ure 7 shows the correlation between Rblr and Paschen 
line luminosities, and Equations (4) and (5) are the best- 
fit results. 
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i?BLR = (50.5±1.0)( 



1042 erg s- 



(0.40±0.01) 

-) It - days, 

(4) 



i?BLR = (50.6 ±1.0) ( 



1042 erg s-i ' 



(0.48±0.01) 



It — days. 

Figure 7 and the above result indicate that the himi- 
nosity of the Paschen hues is closely connected to the 
radius of BLR. 

These tight correlations support the idea that the 
Paschen line emissions originate from the BLR region 
similar to the area where the broad lines of the Balmer 
lines arise. Thus, it should be possible to derive Mbh 
from Pa and P/3 just like Mbh can be estimated from 
H/3 and Ha. Now, we move on to the derivation of Mbh 
estimators based on Paschen lines. 

3.2. Mass Estimators 

In this section, we derive new Mbh estimators 
based on Pa and P/3 lines over the mass range of 
10^-4 ^ 10^-^ Mq. These Mbh estimators are anal- 
ogous to the Mbh estimators based on Ha and H/3 lines 
(Greene & Ho 2005), with Lpa and Lp^ serving as a mea- 
sure of the BLR size, and FWHMs of the Paschen lines 
as the velocity term in the virial mass estimator. Mathe- 
matically, we need find three unknown parameters, a, b, 
and c in the following equation. 



fit with rms scatter of 0.22, 0.23 dex for both Mbh (Pet) 
and Mbh (P/?)- 

The second and the third methods utilize a direct fit 
of the Mbh values. In the second method, we fix the 
exponent of the velocity term c to 2 as expected in the 
virial relation, while setting the overall normalization, a, 
and the exponent of the luminosity term, b, as free pa- 
rameters. For this, we performed a linear regression fit in 
logarithmic scale as in Equation (6), using the REGRESS 
procedure in IDL. 

For the second method, we obtain the following rela- 
tion: 



_ iq7.16±0.04/ ^Pa 

Mq ^ 1042 erg s-i 



M 
M^ 



10 



7.13±0.02 



L 



PI3 



0.49±0.06/ FWHMpa ^ 

403 km s-^' 
(9) 

0.48±0.03, FWHMp^ ,2. 



1042 erg s-i • 



10'^ km s 



(10) 

In Figure 9, we compare the input Mbh versus Mbh 
from the Paschen estimators of the 2nd method. We find 
that the estimators in Equation (9) and (10) reproduce 
Mbh with rms scatters of about 0.20-0.24 dex, not much 
from the method 1. 

In the third method, we treat all the coefficients in 
Equation (6) as free parameters. To derive the coeffi- 
cients in the relation, we use a multiple variable linear 
regression fit with the REGRESS fitting procedure in 
IDL. 

For the third method, we obtain the following relation: 



log(M) = a + b log(i) + c log(FWHM) 



(6) 



The expected values are c = 2 from the virial theorem, 
and b = 0.5 if L ~ ^blr (Dibai 1977). Since there are 
several ways to derive A^bh estimators and the results 
may not be identical, we derive Mbh estimators in three 
different ways as described below. 

First, we derive Mbh estimators by simply replac- 
ing the line luminosities and FWH Ms of the Hydrogen 
Balmer line mass estimators (Al) of (jGreene "fcllol[2007l) 
with those of the Paschen lines using the relation pre- 
sented in Section 3.1. We used the Ha-based Mbh esti- 
mator since the comparison of Ha and the Paschen line 
properties show smaller scatter values than H/3. We re- 
placed the Ha luminosity and FWHM in the base es- 
timator with those of Pa or P/3 lines. A factor of 1.8 
is multiplied to the Ha-based Mbh estimator of Greene 
& Ho (2005) so that the new relation is normalized to 
the virial factor of / = 5.5. The Mbh values for the 
left-hand side of Equation (6) are those listed in Table 
[TJ The Paschen line based Mbh estimators derived this 
way are given below. 



' 10-^ km s^i ' 



(7) 



— = in7.31(- ^Pa 

Mq ^ 1042 erg s-i 



0.48±0.03<- FWHMpa n1.68±0.12 

403 km s-1'^ 



(11) 



M 

Mo 



10"-40( 



1042 erg s- 



^0.46±0.02j^_PWHMp/3 ^i.4i±o.09 



10^ km s 



^ _ 1 r)7.29±0.10/ ^Pa \0.43±0.03(- FWHMpa ^i.92±o. 

Mq~ 4o42 erg s-1^ 



^ = iq7.33±0.10/ Lpp n0.45±0.03(- FWHMp/j .i.69±o 

Mq ^ 1042 erg s-i ^ ^ 103 km s"! ' 

Figure 8 shows Mbh from this estimator versus the 
input Mbh values. The estimators provide a reasonable 



(12) 

Figure 10 shows the comparison of the input Mbh ver- 
sus Mbh from the estimators in Equations (11) and (12). 
The rms scatter in the derived Mbh is reduced to 0.19, 
0.21 dex while the exponents of the velocity term are 
found to be about 1.4-1.7 rather than the value of 2 ex- 
pected in a virial relation. Compared to the 2nd method, 
the third method improves the fitting accuracy by a small 
amount in terms of reducing the rms scatter by 0.03 dex. 

Also, note that the scatter in the Mbh estimator is 
intrinsic (e.g., the uncertainty in the input Mbh val- 
ues) rather than dominated by the measurement errors 
in the FWHM and line luminosity measurements. The 
estimated l-cr errors in FWHM and line luminosity mea- 
surements are 0.04 dex and 0.05 dex respectively, and 
such errors can produce a scatter in the Mbh estimators 
at the level of 0.09 dex only. 

L8 

4. DISCUSSION 

4.1. Consistency between three Mbh estimators 

16 We derived the Mbh estimators in three different ways. 
Ideally, each estimator should give a Mbh value consis- 
tent with the others within the intrinsic rms scatter for a 
given set of realistic FWHM and line luminosity. To see 
if there is any noticeable deviation among the derived BH 
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mass values from the three different estimators, we com- 
pare Mbh values derived from the three estimators. For 
the input sets of FWHM and line luminosities, we convert 
the available FWHM and L(5100) of SDSS quasars 
from Shen et al. (2008) to the Paschen line quantities, 
which gives realistic sets of FWHM and line luminosity 
values in a large range of FWHM and luminosity param- 
eter space occupied by quasars. To convert FWHMh^ to 
FWHMp„ or FWHMp^ and L(5100) to ipa and Lpp, 
we use the the empirical relations presented in Section 
3.1 of this paper. The derived Paschen line FWHM and 
luminosities are fed into the new A/bh estimators. 

Figure 12 shows the comparison of Mbh derived from 
three different estimators. The left panel of the figure 
shows the comparison between Mbh values from the 
method 1 against those from the method 2. There is 
a small systematic offset of 0.1, 0.15 dex in the derived 
Mbh values at Mbh 10^° ° Mq, Mbh 10^ ° M© 
with the overall rms scatter of 0.06 dex. Considering 
that the rms scatter in these two estimators is about 
0.2 dex, these two estimators produce Mbh values well 
within the rms scatter in their relation. The right panel 
of the Figure 12 compares A/bh from the method 2 ver- 
sus the method 3. There is a significant systematic offset 
of 0.2 dex at the high mass end of 10^" Mq and the low 
mass end of 10^ Mq - both are comparable the rms scat- 
ter of the Mbh estimators. The systematic offset results 
mainly from the difference in the exponent of FWHM in 
the A^BH estimators. The derived Mbh is proportional to 
FWHM^ in the method 2, while A^bh from the method 
3 is proportional to FWHM^-^. Therefore, the method 2 
is bound to produce Mbh values larger than the method 
3 at large Mbh while the opposite trend appears at low 
Mbh. 

It is difficult to determine which method gives the best 
estimate of A/bh- The method 3 gives the smallest scat- 
ter, but the FWHM exponent of 1.5 is less physical than 
the index of 2 in the method 2. The estimator from 
the method 3 may be proved to be the best choice, if 
a good explanation can be provided why the exponent 
of FWHM is 1.5, not 2. The Mbh estimator from the 
method 3 has a resemblance to the fundamental plane 
relation of early-type galaxies (Dressier et al. 1987; Djor- 
govski 1987) where it still remains controversial why the 
relation is slightly tilted off from the virial relation (e.g., 
see Jun & Im 2008), so a similar physical mechanism 
that produces the tilt in the fundamental plane might be 
in work for the Mbh estimator. For now, we prefer the 
method 2 the most simply because the formula can be 
justifiable easily on a physical basis. 



The line ratios of H/3 through Pa lines can give us clues 
on the physical conditions of the BLR region. The line 
ratios can also serve as a basis for determining the ex- 
tinctions at different wavelengths in BLR of dusty AGN. 
Therefore, we present the line ratios of the Hydrogen 
lines of our sample. 

Figure 11 shows the line ratios of H/3, Ha, P/3 to 
Pa of our sample. Overplotted are theoretical line ra- 
tios based on the co mputation from the CLOUDY code 
(jFerland et al.l[T998t ). Using the CLOUDY code, we cal- 
culate the expected line ratios by varying three parame- 
ters, the shape of the ionizing continuum (a=-1.0, -1.5), 
the ionization parameter {U = 10°'^, 10~°-^, 10"-^-^) 
and the hydrogen density {n — 10^, 10^", 10^^ cm^'^). 
We find that the observed line ratios are reproduced 
most successfully with a set of parameters, a = —1.0, 
U = 10~^'^ and n = 10^ cm"'^. Some line ratios that are 
below the theoretical expectation can be explained nat- 
urally if there is a small amount of extinction with color 
excess of E{B — V)< 0.4 mag. The mean line ratios are 
found to be = 2.70, ^ = 9.00, and ^ = 0.91. These 
values can be used as a basis for determining extinctions 
in the BLR of dusty AGN. The perfect fit of the line ra- 
tios as a function of wavelengths is probably not possible 
with a simple model, considering that the broad emission 
lines at different wavelengths are not likely to originate 
from exactly the same BLR regions as we discussed in 
Section 3.1. Nevertheless, the observed line ratios are 
consistent with the expected ratios from a single set of 
ionization parameters. 

5. CONCLUSION 

We derived new Mbh estimators, using Hydrogen Pa 
and P/3 lines of Type-1 AGNs. The derived estimator al- 
lows the determination of A/bh at the accuracy of ~ 0.2 
dex, and they will be useful for estimating Mbh of red, 
dusty AGNs. Our analysis of the Paschen lines with re- 
spect to Ha and H/3 lines shows that the luminosities and 
FWHMs of the broad components of the Paschen lines 
correlate well with those of the Balmer lines. The Hydro- 
gen line ratios from H/J through Pa are consistent with a 
Case B recombination with the parameters of a = —1.0, 
U = 10~^'^ and n = 10^ cm~'^. The mean line ratios are 
= 2.70, ^ = 9.00, and ^ = 0.91 which can be used 

Pa . ' Pa ' Pa 

to estimate for the amount of dust extinction present in 
red AGNs in future. Future applications of these results 
on red, dusty AGNs will enable us to better understand 
the nature of red, dusty AGNs whose properties are still 
hidden behind a wall of dusty gas. 



4.2. Line Ratios 
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TABLE 1 

Paschen Parameters And BH Masses 



AGNs Name 


Redshift 


BHmass 






i'f'' 








i'o 






z 


log(M/M0) 


FWHM 


AFWHM 


Flux 


Aflux 


1 WHM 


A "m T T'T T TV T 

AFWHM 


Flux 


Aflux 






n — li 
[km s J 


n — li 
[km s J 


r —1 — 2i 
[erg s cm J 


r —1 — 2i 
[erg s cm J 


n — li 
[km s J 


n —11 
[km s J 


r —1 — 2i 
[erg s cm J 


r —1 
[erg s cm 


30273=' 


0.158 


8.94': 


2608 


— 


6.95E-13 


— 


2635 


— 


8.33E-13 


— 


Mrk876'' 


0.129 


8.44': 


5433 


- 


7.81E-14 


- 


4965 


- 


1.31E-13 


- 


PG0844+349'' 


0.064 


7.96=: 


2144 


- 


1.06E-13 


- 


1967 


- 


1.21E-13 


- 


MrkllO* 


0.035 


7.40'= 


1689 


- 


1.50E-13 


- 


1604 


- 


1.99E-13 


- 


Mrk509*' 


0.034 


8.15"= 


2755 




5.20E-13 




2747 




6.60E-13 






U.Uoo 


ft 1 7C 


4D1U 




o.oiytlj-lo 




40 (54 




a OflTT 1 "i 
D.yoCj-iO 




Mrk817*' 


0.031 


7.59c 


4988 




1.39E-13 




4202 




1.39E-13 




Mrk335'' 


0.026 


7.15'= 


1633 




1.67E-13 




1431 




1.40E-13 




Mrk590'' 


0.026 


7.67'= 


3565 


_ 


1.46E-14 


_ 


4259 


_ 


2.82E-14 


_ 


NGC5548'' 


0.017 


7.82'= 


5891 


- 


1.77E-13 


- 


4102 


- 


1.64E-13 


- 


NGC4151'' 


0.003 


7.12'= 


4204 


- 


1.21E-12 


- 


- 


- 


- 


- 


H1821+643'' 


0.297 


9.44'' 


4717 


- 


1.92E-13 


- 


- 


- 


- 


- 


PDS456'' 


0.184 


8.70^ 


1828 


- 


3.79E-13 


- 


1785 


- 


4.70E-13 


- 


Mrk290» 


0.030 


8.09't 


3818 


- 


1.55E-13 


- 


3180 


- 


1.64E-13 


- 


H2106-099'' 


0.027 


7.72'^ 


2155 


- 


1.70E-13 


- 


1550 


- 


2.35E-13 


- 


Mrk79'' 


0.022 


7.94'^ 


3163 


- 


1.86E-13 


- 


3051 


- 


2.58E-13 


- 


NGC4593^ 


0.009 


7.83'^ 


3407 


- 


2.87E-13 


- 


- 


- 


- 


- 


NGC3227^ 


0.004 


7.42'^ 


2680 


- 


3.48E-13 


- 


- 


- 


- 


- 


HE 1228+013" 


0.117 


7.87't 


1731 


- 


6.62E-14 


- 


1725 


- 


7.46E-14 


- 


H1934-063'' 


0.011 


6.80'' 


1241 


- 


1.37E-13 


- 


1213 


- 


1.77E-13 


- 


IRAS1750+508'' 


0.300 


8.33'' 


1758 


- 


4.50E-14 


- 


- 


- 


- 


- 


SDSSJ000943. 1-090839. 2b 


0.210 


8.44'= 


4817 


818 


8.41E-15 


1.84E-15 


3841 


724 


8.31E-15 


2.10E-15 


SDSSJ005812. 8+160201. 3*^ 


0.211 


8.29'= 


3158 


409 


1.73E-14 


1.96E-15 


2775 


357 


2.23E-14 


2.37E-15 


SDSSJ010226. 3-003904. 6*^ 


0.295 


7.67'= 


1565 


184 


2.42E-14 


1.87E-15 


1234 


142 


2.86E-14 


2.04E-15 


SDSSJOllllO. 0-101631. 8*^ 


0.179 


8.20'= 


3695 


496 


1.48E-14 


1.77E-15 


2617 


380 


1.41E-14 


1.99E-15 


SDSSJ015530. 0-085704.0'' 


0.165 


8.85'= 


7355 


1032 


1.53E-14 


2.38E-15 


5074 


737 


1.36E-14 


2.24E-15 


SDSSJ015910. 0+010514.5'' 


0.217 


7.98'= 


2657 


404 


l.OlE-14 


1.48E-15 


2574 


455 


1.02E-14 


1.89E-15 


SDSSJ015950. 2+002340. 8*' 


0.163 


8.09'= 


2657 


391 


2.65E-14 


2.53E-15 


2716 


627 


3.78E-14 


3.99E-15 


SDSSJ021707.8-084743.4'' 


0.292 


7.84'= 


1970 


552 


3.49E-15 


1.14E-15 










SDSSJ024250. 8-075914. 2^ 


0.378 


8.00'= 


1329 


194 


9.28E-15 


1.15E-15 










SDSSJ031209.2-081013.8'' 


0.265 


8.29*= 


3540 


620 


1.12E-14 


2.01E-15 










SDSSJ032213.8+005513.4'' 


0.185 


8.03" 


2014 


234 


2.58E-14 


1.89E-15 


1894 


220 


3.40E-14 


2.47E-15 


SDSSJ150610.5+021649.9'' 


0.135 


8.11" 










2632 


330 


2.26E-14 


2.54E-15 


SDSSJ172711.8+632241.8'' 


0.218 


8.61" 










5087 


1192 


1.13E-14 


2.86E-15 


SDSSJ211843.2-063618.0t' 


0.328 


8.42" 


2485 


430 


6.58E-15 


1.29E-15 










SDSSJ234932.7-003645.8'' 


0.279 


8.19" 


2578 


360 


1.02E-14 


1.47E-15 










SDSSJ235156.1-010913.3'' 


0.174 


8.89" 


6137 


705 


4.42E-14 


3.53E-15 


4946 


553 


5.14E-14 


3.65E-15 



^ FWHM and Flux of these objects come from L08. 

^ FWHM and Flux of these objects measured from spectra presented in GOG. 

^ The Mbh values arc determined from reverberation mapping techniques (Vestergaard & Peterson 2006). 

The Mbh values are based on H/3 widths and luminosities in LOS. 
^ The Mbh values are based on H/3 widths and L(5100) derived from our fit to SDSS spectra. 
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TABLE 2 

Coefficients of the Empirical Correlation between Balmer and Paschen Lines 



Num 


Y 


X 




FWHM 






Line luminosity 




A 


B 


rms (dex) 


C 


D 


rms (dex) 


1 


Ha 


Pa 


0.074±0.038 


0.934±0.084 


0.045 


0.916±0.014 


0.961±0.025 


0.141 


2 


H/3 


Pa 


0.105±0.037 


1.017±0.080 


0.057 


0.444±0.013 


0.910±0.018 


0.188 


3 


Ha 


P/3 


0.076±0.038 


0.821±0.075 


0.051 


0.985±0.012 


L008±0.015 


0.117 


4 


H/3 


P/9 


0.113±0.033 


0.895±0.068 


0.058 


0.517±0.011 


0.943±0.013 


0.162 
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Wavelength (/xm) Wavelength {/J,m) 



Fig. 1. — An example of the continuum fitting of the GOG NIR spectra. On the top panel, we show the NIR spectrum of 
SDSSJ032213. 8+005513. 4, where the blue solid line indicates the continuum fit to the Po wavelength region and the red solid line shows 
the continuum fit to the P/3 wavelength area. The wavelength ranges that were used for determining the continuum are indicated with the 
vertical dashed lines. In the bottom panels, we show the expanded views around P/3 and Pa lines together with the determined continuum. 
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Fig. 2. — Results of the fitting of Pasctien lines with a Gaussian function to the GOG sample. AGNs with both Pa and P/3 measurements 
are shown here. The ordinate is in units of 10~^^ erg s~^ cm"-^ fj,m~^. The continuum is already subtracted. The black solid line indicates 
the best-fit and the grey solid line indicate observed spectrum. 




Fig. 3. — Gaussian fits of the Pasclien lines of the GOG sample. AGNs with only Pa or P/3 measurements are shown here. The meaning 
of the lines is identical to Figure 3. 
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Fig. 4. — The comparison of the P/3 line profiles against the Hq line profiles. Only objects with P/3 line S/N > 25 are shown here 
(histogram with errors). The Ha lines are indicated with thick solid lines, and the thin dotted lines and the thick dashed lines represent 
the broad and the narrow components of the Ho lines. The P/3 profiles are normalized to have the maximum value of 1 while the Ho are 
normalized to have the same total fiux as the normalized P/3 profiles. The figure shows that the P/3 line profile shapes are consistent with 
the Ho profile shapes. 
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Fig. 5. — The comparison of the FWHM widths of the Paschen Hnes versus the Balmer Hnes. The top panels compare the Po FWHMs 
versus the Balmer line FWHMs, while the bottom panels show the comparison of the P/3 FWHMs versus the Balmer line FWHMs. The 
open circles are for AGNs from G06, while the red triangles are for the L08 sample. The solid line indicates a line where the Balmer and 
the Paschen quantities are identical. The dotted line indicates the best-fit line between the two quantities as described in the text. 
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Fig. 6. — The comparison of the Paschen hne luminosities against those of the Balmer Hnes. The top two panels compare the Pa fluxes 
against the Balmer lines fluxes, while the bottom two panels compare the P/3 fluxes against the Balmer line fluxes. The symbols are same 
as those in Figure 5. The meaning of dotted line is identical to that in Figure 5. 
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Fig. 7. — The correlation between the Pa (left) and the P/3 (right) luminosities and the BRL sizes. The BLR sizes are derived from the 
5100A luminosities. The good correlation (albeit over 40 It-days to 110 It-days) suggests that the Paschen line luminosities can be used to 
estimate the BLR size. The solid lines indicate the best-fit relation of the correlation. 
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Fig. 8. — Comparison of the Paschen-line based Mbh versus Mbh derived from the reverberation mapping or optical spectra using the 
1st method described in the text. Black circles indicate the GOG sample. Blue crosses show the LOS sample with Mbh derived from the 
reverberation mapping ( Vestergaard & Peterson 2006). Red squares are for the LOS sample with Mbh estimated from single epoch optical 
spectra. The solid line indicates an one to one function with Mbh based on existing method and it upon properties of Paschen line. 
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Fig. 9. — Similar to Figure 8, but with the Paschen-based Mbh where Mbh are derived with an estimator for which the power-law index 
of FWHM is fixed to 2. The meaning of the symbols and line are identical to Figure 8. 
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Fig. 10. — Similar to Figure 8, but with the Paschen-based A/bh where Mbh are derived with an estimator for which the power-law 
index of FWHM is set free during the fit. The meaning of the symbols and line are identical to Figure 8. 
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Fig. 11. — Line ratio of type-1 AGNs compared with theoretical models from CLOUDY code. 
These models show that typical type 1 AGN condition is a= -1.0, n= lO" cm""^ and U= 10~^'^. 
from various condition, within a= -1.5 ~ -1.0, n= lO' ~ 10^^ cm^'^ and U= 10~^'^ ~ lO'''^. 



Fluxes are normalized to that of Pa. 
The grey lines are theoretical models 
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Fig. 12. — The left panel compares Mbh from the method 1 and the method 2. In the right panel, we compare Mbh fr'om the method 
2 versus the method 3. The Mbh values agree well with each other in this case. The data points come from SDSS DR5 quasars in Shen 
et al. (2008) where the Paschen FWHM and luminosities are estimated from the corresponding H/3 quantities. The solid line indicates the 
case where the BH masses derived from different methods are identical. 



